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Angle-resolved photoemission spectroscopy (ARPES) is used to study the spectral function of the 
optimally doped high-Tc superconductor (Bi,Pb)2Sr2CaCu20g+i in the vicinity of the antinodal 
point in the superconducting state. Using a parameterized self-energy function, it was possible to 
describe both the coherent and the incoherent spectral weight of the bonding and the antibonding 
band. The renormalization effects can be assigned to a very strong coupling to the magnetic res- 
onance mode and at higher energies to a bandwidth renormalization by a factor of two, probably 
caused by a coupling to a continuum. The present reevaluation of the ARPES data allows to come 
to a more reliable determination of the value of the coupling strength of the charge carriers to the 
mode. The experimental results for the dressing of the charge carriers are compared to theoretical 
models. 
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I. INTRODUCTION 

The dressing of the charge carriers in high-Tc super- 
conductors (HTSCs) is still one of the most exciting top- 
ics in solid state physics. The HTSCs are a paradigm 
for the transition of a correlated system from an insu- 
lating to a metaUic state. The dressing of the charge 
carriers in HTSCs is Hkely caused by the same interac- 
tion that gives rise to the superconductivity, hence the 
understanding of the quasiparticle self-energy may help 
to understand the origin of the mechanism of high-Tc su- 
perconductivity. The dressing can be studied by various 
experimental methods but angle-resolved photoemission 
spectroscopy is the only method which gives a quantita- 
tive information on the momentum dependence of those 
renormalization effects. In HTSCs there are two im- 
portant regions on the Fermi surface: the nodal region, 
where the diagonal of the Brillouin zone cuts the Fermi 
surface and where the d— wave superconducting order pa- 
rameter changes sign. This region mostly contributes to 
the transport properties, particularly in the underdoped 
region, where a pseudogap opens up, squeezing the Fermi 
surface to a region near the nodes. The other (antinodal) 
region is one where the edge of the Brillouin zone cuts the 
Fermi surface. In this region, the order parameter in hole 
doped superconductors has its maximum. This region is 
therefore mostly relevant for the studies of the supercon- 
ducting properties. There are numerous ARPES studies 
on the renormalization effects near the nodal point (i*2ii 
but only a few studies are concentrated at the antinodal 
point '^i^i^i'^i^i'^. 



In the bilayer systems the study of the antinodal point 
is complicated by the bilayer splitting, which could not 
be resolved for 15 years. On the other hand, only in 
the bilayer system of the Bi-HTSC family the entire su- 
perconducting region from underdoped (UD) via opti- 
mally doped (OP) to overdoped (OD) can be studied. 
In the superconducting state a well pronounced peak- 
dip-hump structure has been detected^. This struc- 
ture was originally explained^iift solely in terms of a cou- 
pling to a bosonic mode, similar to the McMillan-Rowell 
explanation of the tunnelling spectra in conventional 
superconductorsii. Later on, it was established that this 
peak-dip-hump structure is partially caused by the bi- 
layer splittingiSiil. By varying the photon energy hv 
in the ARPES experiments, and exploiting the different 
energy dependence of the matrix elements for the excita- 
tions from the bonding and the antibonding bands, it be- 
came possible to separate the two bandsStiSiiSiiiii^ and to 
extract the full energy- and momentum-dependent spec- 
tral weight separately in each of the ba nds. This proce- 
dure allowed the authors of R,efs. ll2ll^ to find the intrin- 
sic peak-dip-hump structure, and to demonstrate that 
the strength of this intrinsic effect is doping-dependent, 
and decreases in going from UD to OD materials. 

An important characteristic of the interaction be- 
tween fermionic and bosonic excitations is the energy- 
dependent, dimensionless coupHng Xe. In theories where 
the fermionic self-energy depends on energy, E, much 
stronger than on the momentum k — kp, this dimension- 
less coupHng is related to the self-energy via T,{E) = 
—EXe- It is also relevant whether the measurements 



are performed in the normal or in the superconducting 
state. We will label the corresponding couplings as A„^_e 
and Xsc,E, respectively. 

If the normal state is a Fermi liquid, Xn,E=o ~ Ki 
is finite, and is often called a dimensionless coupHng 
constant. It determines the mass renormalization of 
the fermionic quasiparticles via m* = m(l + A„). The 
coupling constant can, in principle, be extracted from 
ARPES measurements of the quasiparticle dispersion in 
the normal state at the lowest energies, however this pro- 
cedure requires one to know both kp and the bare mass, 
m. In previous analysis^, the mass, in, was extracted 
from a tight-binding model with parameters derived from 
a fit of the Fermi surface and from the quasiparticle dis- 
persion measured along the nodal directioniSiii. The 
analysis of the experimental data in the antinodal region 
yielded A„ ^ 1.5 both in UD and OD materials. This re- 
sult should be contrasted with values^^ of A„ < 1 at the 
nodal point. It is consistent with expectations as for non 
rotationally-invariant systems the coupling A,i depends 
on the position on the Fermi surface. 

In the superconducting state, the measured quasipar- 
ticle energy in the antinodal region is bounded by the su- 
perconducting gap. A, and it becomes an issue at which 
energy one extracts the coupling Xsc,e from the data. In 
previous analysis, the coupling was extracted from the 
self-energy measured at \E\ > A. This coupHng Xsc,a 
turns out to be larger than A„, and it also rapidly in- 
creases from OD to UD samples (Asc,a ~ 8 for dopant 
concentration 0.12). 

In this communication we extend our previous analy- 
sis of the antinodal self-energy in the superconducting 
stated and show how one can extract the coupling at 
zero frequency Xsc,e=o = Xsc from the ARPES data. 
We find that Xsc is smaller than Asc,A and within a cer- 
tain model is also smaller than the normal state cou- 
pling A„, in agreement with earlier calculationsi^. We 
show that the large value of Xsc, a and its strong dop- 
ing dependence are at least partially due to the fact that 
the fermionic self-energy in a superconductor actually di- 
verges at \E\ = A -I- rio, where flo is the energy of the 
bosonic mode. If the bosonic mode is the spin resonance 
peak, its energy decreases with decreasing doping. Then 
\E\ = A and \E\ = A + flo come closer to each other 
in the UD regime, and Xsc, a strongly increases. This is 
consistent with the analysis in Ref. H. 

Our present analysis is based on the measurements of 
the quasiparticle spectral function in the antinodal region 
of the high-Tc superconductor (Bi,Pb)2Sr2CaCu208+5 
(BiPb2212) in the superconducting state. We go be- 
yond a previous ARPES study which has analyzed the 
energy dependence of the spectral weight just at the (tt, 0) 
poin1j2fi, and study the whole antinodal region. We inter- 
pret our results in the superconducting state in terms 
of model self-energy function which is composed of two 
terms. The first and dominant term is due to a strong 
coupHng of the charge carriers to a single bosonic mode. 
The second term describes a band renormalization at 



higher energies and is assumed to have a Fermi-liquid 
form. We extract both couplings from the fits to the 
data. We used two models for electron-boson coupling. 
The first model is a one-mode model for an interaction 
with an Einstein boson, which is assumed to be inde- 
pendent on fermions. Second is a collective mode model, 
in which the bosonic spectrum in the normal state is 
rather fiat and incoherent, but splits into a mode and into 
gapped continuum in the superconducting state due to 
the feedback effect from the pairing. This second model 
is appropriate if the boson is a spin collective mode of 
fermions. We obtain a rather good agreement between 
the parameters derived from the analysis of the experi- 
mental data using the model self-energy function and the 
calculated values using the collective mode model. This 
yields a strong indication that the dominant part of the 
renormalization of the fermionic dispersion is due to a 
coupHng of coHective spin excitations. 

The paper is organized as follows. In Sec. II we re- 
view the two fermion-boson models in the normal and 
the superconducting state. The experimental setup is 
discussed in Sec. III. In Sec. IV we present the exper- 
imental results together with the data analysis. In Sec. 
V we discuss the results and compare them with other 
renormalization effects studied by ARPES in solid state 
physics. The conclusions of our study are presented in 
Sec. VI. 



II. 



THE FERMION-BOSON MODELS 



The coupling of the charge carriers to bosonic exci- 
tations is the minimum model to understand the spec- 
tral function of the HTSCs at the antinodal point. We 
start with an assumption that the Fermi energy Ep is 
much larger than the mode energy flo. The validity of 
this assumption for very underdoped cuprates has been 
questioned recentl y^^'^^ because there the bandwidth is 
strongly reduced due to correlation effects associated 
with Mott physics. Here we restrict our analysis to near- 
optimally doped cuprates for which there is little doubt 
that Ep » r^o since Ep ~ 1 eV in this case. 

Both fermion-boson models have been discussed ear- 
lier in the literatureiSiSSiSSiSiiSSiSS. We review them here 
again in order to specify the parameters which can be de- 
rived from ARPES. We also present several new results 
for the collective excitations model. 

The dynamics of an electron in an interacting system 
can be described by a Green's functional 



G{E,k) 



1 



E-ek~^{E,ky 



(1) 



where Y.{E,k)=T,'{E,k)+iT,"{E,k) is the complex self- 
energy function which contains the information on the 
fermion-boson interaction, and Ck is the bare quasiparti- 
cle dispersion. Near the Fermi surface e^ = vp{k — kp), 
where vp = kp/m, and m is the bare mass. It is custom- 
ary to use the tight-binding form for efc. 



ARPES experiments measure the product of the spec- 
tral function A{E,k), the Fermi function, and a tran- 
sition matrix element, convoluted with the experimental 
resolution. The spectral function is related to the Green's 
function a o^^i^^ 



A{E,k) 



-ImG{E,k) 
1 



£"(£;, /c) 



iT[E-ek-^'{E,k)Y + [i:"{E,k)\ 



■(2) 



For S = 0, i. e., for the non- interacting case, the spectral 
function A{E, k) = 5{E - efe). 

For the description of the spectral function in the su- 
perconducting case, two excitations have to be taken into 
account: the electron-hole and the pair excitations. This 
transforms the Green's function into a (2x2) matrixj^ii or, 
equivalently, to the emergence of normal and anomalous 
components of the Green's function. Accordingly, the 
self-energy also has a normal part S(i?, k) and anoma- 
lous part $(£', k). The two self-energies are related to E', 
the renormalization function Z(E,k), and to the super- 
conducting gap A(_E, k) via 

E - E{E, k) = EZ{E, k), ^{E, k) = Z{E, k)A{E, k) 

(3) 
In general, the superconducting gap A(i?, k) is also a 
complex function and depends on both parameters. The 
energy dependence is not crucial, though^!, and we just 
neglect it for simpHcity, i.e., replace a complex A{E,k) 
by a real A(A:). The momentum dependence of A(fc) is 
that of a dx2_y2 gap. In the antinodal region, the gap 
is near its maximum, its momentum dependence is weak 
and we will neglect it as well, i.e., further approximate 
A(fc) by A and Z{E,k) by Z{E). The spectral function 
is then given bji24 

Using our definition of the coupling constant, 

ZiE) = l + X,c,E. (5) 

Below we consider two models for electron-boson in- 
teraction. In the one-mode model we define the self- 
energy due to the coupling to a single bosonic mode 
as S' = — A^S with Ag = A**. In the collective mode 
model, we treat the renormalization due to a distribution 
of bosonic modes. The corresponding coupling constant 
is called A"^. 



A. One-mode model 

In the one-mode model, it is assumed that electrons 
interact with an Einstein boson whose energy is Oq in- 
dependent on whether the system is in the normal or 
in the superconducting state. For the normal state the 
mechanism leading to a finite lifetime of a photohole is 



Table I: Bosonic excitations which couple to the charge carri- 
ers together with their characteristic energies in HTSCs 



system 



excitations characteristic energy(meV) 



ion lattice phonons 

spin lattice/liquid magnons 
e-liquid plasmons 



90 

180 

1000 





Figure 1: Bosonic scattering mechanism which contributes to 
the imaginary part of the self-energy, (a) normal state; (b) 
superconducting state. 



illustrated in Fig. QJa). The hole is filled by a transition 
from a state at lower binding energy via an emission of a 
bosonic mode. Such bosonic excitations may be electron- 
hole excitations, phonons, spin excitations, plasmons, ex- 
citons etc. Relevant excitations for HTSCs are listed in 
Table I together with their characteristic energies. 

For a constant density of states and for the tempera- 
ture T = 0, the fermionic self-energy is given bySi 



^{E) = ^ J dE'x{E') J dtuG{E + E', k) 
where x(-E') is the bosonic propagator 



X{E') 



no 



fil - E'-^ - iS 



(6) 



(7) 



In the normal state, J dekG{E + E',k) = — i7rsgn(£' -|- 
E'). Substituting this into lO and separating real and 
imaginary parts of the integral, we find, for E < 



^"iE) = ^x'„noe{\E\^no). 



(8) 



T."{E) is zero up to the absolute value of the mode energy 
Oq . This is also clear from Fig.Q]since the photohole can 
only be filled when its binding energy is larger than fig. 
At \E\ > Qo, ^"{E) is a constant (see Fig.[l(b)). S'(£;) 
shows a logarithmic singularity at the mode energy, Qq 
(see Fig. [21(a)). At low energies there is a linear energy 
dependence of S' and the negative slope —d'S'/dE = Afj 
determines the coupling constant at zero energy. 

In Fig. [21 (c) and (d) we have plotted the renormal- 
ization function Z(E) for the same parameters. The real 
part shows again a singularity at fig and a constant value 
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Figure 2: Real (a) and imaginary (b) part of the self-energy 
function and real (c) and imaginary (d) part of the renormal- 
ization function for a coupling to a mode at Q.q = 40 meV and 
a coupling constant X„=8. Solid line: normal state, dashed 
line: superconducting state with a gap A = 30 meV 



at zero energy. This value minus 1 again determines the 
coupHng constant A^. Within this model the same cou- 
pling constant can also be obtained from the measure- 
ments of the quasiparticle linewidth at large negative en- 
ergies as X'^ = —21] (— oo)/(7rno)- S (— oo) is the step 
height of S (E) aX E = ^Iq in the one-mode model. Hence 
both A^ and S (— oo) are the measures of the coupling 
strength to the bosonic mode. Consequently, together 
with the mode energy Oq, both can be used to determine 
the self-energy in the one-mode model. 

In Fig. El (a) and (b) we have displayed the calculated 
spectral function in the one-mode model for Afj = 1 and 
A^ = 8, respectively. Compared to the bare particle dis- 
persion, efc, given by the red dashed line, for |E| < fla 
there is a mass renormalization, i.e., a reduced dispersion 
and no broadening, except the energy and momentum 
resolution broadening, which was taken to be 5 meV and 
0.005 A~^, respectively . For |E| > flo, there is a disper- 
sion back to the bare particle energy. Moreover, there is a 
broadening due to a finite S", increasing with increasing 
A^. For large A^, the width for constant E scans is, at 
least up to some energy, larger than the binding energy 
of the charge carriers and therefore they can be called 
incoherent in contrast to energies |E| < flo or very high 
binding energies, where the width is smaller than the 
binding energy and therefore the states are coherent. 
The change in the dispersion is often termed a "kink" 
but looking closer at the spectral function, in particular 
for high A^ , there is a branching into two dispersion arms 
touching each other at the branching energy ¥13=^0- 
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Figure 3: Calculated spectral function A(E,k) for a coupling 
of the charge carrier to a mode with an energy no=40 meV 
(a) and (b) normal state, (c) and (d) superconducting state 
with a superconducting gap, A=30 meV. (a) and (c) coupling 
constant A,i = l, (b) and (d) coupling constant A„=8 



The following information can be obtained from a one- 
mode model spectral function A{E, k) in the normal 
state. When performing constant-E scans, very often 
called momentum distribution curves (MDCs), one ob- 
tains Lorentzians. The maximum observed in the MDCs 
determines the renormalized dispersion. By comparing 
this dispersion close to Ep to the bare particle disper- 
sion one can extract the coupHng constant Aj'j. The 
width of the Lorentzians for \E\ < flo in this model 
should be determined by the energy and momentum res- 
olution. For \E\ > rio the resolution effects are small 
compared to the intrinsic width, W, and one can derive 
S"(— 00) ~ ~vpW/2 where vp is the bare-particle Fermi 
velocity. From the onset of a finite intrinsic width, one 
obtains the mode energy Qq. This parameter also can 
be obtained by constant-E cuts, very often called energy 
distribution curves (EDCs). Looking at Eq. (2), one re- 
alizes that for large e^, i.e., far away from the Fermi wave 
vector, kp, the spectral function around E ^ ilo is deter- 
mined by S". The edge of such a cut determines again 
no. 

In the superconducting state, the self-energy for i? < 
is still given by ©, but the fermionic Green's function 
now has the form 



Gsc{E,k) = 



E + €k 



E^ 



iS 



(9) 



Substituting this into (JSJ , evaluating the integral over e^ 
and separating real and imaginary parts, we obtain for 



E <0: 



one to neglect the energy dependence of XscE at |i?| > A, 






2 " ^(£; + f7o + A)(S + f7o-A) 



(10) 



In Fig.|2lwe plot S(£') and ^(-E) for the superconducting 
state with A = 30 meV. A small 5 has been used to reduce 
the singularities. One realizes that due to the opening of 
the gap the singularities of S' and Z' are shifted to higher 
binding energies and that the edge in E" transforms into 
an overshooting edge. At vanishing E^ ^'(E) is still linear 
in E, but due to the shift of the singularity to higher 
binding energies the slope is now reduced and given by 






0) 



E 



= A 



dx 



1 



(x2 + 1)3/2 l+x2(A/f7o)2 
(11) 

The reduction of the coupling constant in the supercon- 
ducting state is concomitant with a reduction of Z{0) 
since Z{0) = 1 + A. 

At k = kp the photohole can only be excited when 
its binding energy is exactly A, or when ImZ{E, kp) ~ 
—ImJ^{E,kp)/E is non-zero. From Fig.Q](b) it is clear 
that for zero temperature and in the clean limit E" or 
the scattering rate is different from zero only when \E\ > 
A -I- rig. This result is obtained from an evaluation of 
Eq. ltTn)l . This implies that i? = — A is separated from 
the region where E {E, kp) is non-zero and therefore the 
spectral function contains a iS— functional peak at £^ = 
—A, and then it becomes nonzero at i? < — (A -I- Qq). 

In Figs.|Hl(c) and (d) we show for the one-mode model 
the calculated spectral function in the superconducting 
state using the same energy and momentum resolutions 
and the same mode energy as before. The gap was set to 
A = 30 meV. One clearly realizes the BCS-BogoHubov- 
like back-dispersion at the gap energy A and besides this, 
a total shift of the dispersive arms by the gap energy. 
Thus the branching energy E^ occurs at —{flo + A). 

The renormaHzed dispersion is obtained from the po- 
sition of the MDC peak of the spectral function. In the 
normal state, the peak position is where the real part of 
G~^{E,k) vanishes. In the superconducting state, there 
is an extra complication due to the fact that Cfc is present 
both in the denominator and in the numerator of the 
spectral function. Like in an earlier studj>i^ we avoid 
this complication and extract the renormaHzed disper- 
sion from 



Cfe 



-ReZ{E)^/E^-A^ 
ReE(^) 
E 



Ve^ - A2. (12) 



In conventional superconductors, the mode energy is 
much larger than the gap. In this situation, Eq. lfTT|l 
yields A^, = a:; (l + O ((A/r!o)2log A/r!o)) « A^ Fur- 
thermore, the same small parameter A/f]o also allows 



such that A 



sc,E 



Kc ~ A„. In this situation, ReZ{E) 



1 -f Afj , and hence the maximum of the spectral function 
is located at 



E=-JA^ + el/{l + XiY 



(13) 



For HTSCs, the gap is comparable to the mode energy 
and therefore Ea. ljl3|l is no longer valid, and the full Eq. 
II12|I should be used to fit the dispersion. There are two 
key differences with Eq. lfT3JI . First, the zero-energy 
values Ajc and Af^ are different. For A = 30 meV and 
Q.0 = 40 meV, i.e., A/fio = 3/4, we obtain from lfTT|l 
Agg = 0.74AJ',. Second, the energy dependence of A^^ 
becomes relevant. Indeed, by analyzing ifTTijI one finds 
that S {E) is discontinuous at i? = — (A + f^o) and di- 
verges as a square-root at approaching E — — (A + i^o)"^ 
When A and 0.^ are comparable, this divergence affects 
the self-energy already at E k. —A. For the parameters 
that we choose, the effect is not large: evaluating the 
real part of the self-energy at i? = — A from (ll()|l we find 
'^sc.A ~ 1-lAsc. However, the effect increases once Oq 
gets smaller. 

Measuring an EDC at kp with high resolution, one 
would expect a peak at A, followed by a region of near- 
zero spectral weight and a threshold of the incoherent 
spectral weight, which appears at fio + A. Such an 
energy distribution is well known from tunnelling spec- 
troscopy in conventional phonon superconductors, except 
that there f^o is often much larger than A. At deviations 
from kp, the peak disperses to larger frequencies while 
the onset of the incoherent spectral weight remains at 
Q.0 -\- A. Once the peak disperses close to VIq -(- A, only 
the threshold at this energy remains visible. 



B. Collective mode model 

For definiteness, we consider the model with the inter- 
action between fermions and their spin collective excita- 
tions with momenta near Q = (7r,7r). The momentum Q 
connects Fermi surface points within antinodal regions, 
and hence antinodal fermions are mostly involved in the 
scattering of nearly antiferromagnetic spin fiuctuations. 

The physics of electron-boson interaction is somewhat 
different in the one-mode and collective mode scenarios. 
Like we said, in the one-mode formalism, one assumes 
that bosons are propagating excitations with a frequency 
f2o, independent on whether fermions are in the normal 
or in the superconducting state. In the collective mode 
model, bosons are Landau-damped in the normal state, 
and their spectral function is described by a continuum 
rather than by a mode. In the superconducting state, 
the low-energy fermionic states in the antinodal regions 
are gapped, and the continuum of bosonic states with 
momenta near Q appears only above the gap of 2 A. In 
addition, the residual attraction between fermions in a 
dx2-y2 superconductor leads to the development of the 



resonance peak at a frequency flo below 2 A. In the OD 
regime, Qq is only slightly below 2A, and the resonance 
is weak. In the UD regime, the resonance frequency de- 
creases. In bilayer systems, such as Bi2212, there are 
two resonances, in the even and in odd channel. The 
resonance frequency in the even channel should vanish at 
the point where the magnetic correlation length diverges. 
The resonance in the odd channel remains finite at this 
point, and, very likely, transforms into the gapped spin- 
wave mode in the antiferromagnetically ordered state. 

The self-energy within the collective mode model has 
been analyzed in Refr.llflland in earlier publications. Be- 
low we briefiy review the existing results and also present 
several new formulas. For definiteness, we consider the 
case of a flat static susceptibility near Q, i.e., assume 
that in the normal state, the dynamical spin susceptibil- 
ity (the bosonic propagator) can be expressed as 



x{E,q)^x{E)^ 



XQ 



1 — lE /ujgf 



(14) 



where uogf is the typical relaxational frequency of spin 
fluctuations. An advantage of using the flat static spin 
susceptibility is that all computations can be done ana- 
lytically. Similar results are also obtained using Ornstein- 
Zernike form of the static susceptibilitj«2i and in FLEX 
computations for the Hubbard model^^. 

In the normal state, the fermionic self-energy due to 
interaction with the gapless continuum of spin excitations 






-X'^LOsf arctan 
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--A>,/ln(l 
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(15) 



where A^ is the dimensionless coupling constant in the 
normal state for the collective mode model. In dis- 
tinction to the one-mode model, the self-energy in Ijl5|l 
has no threshold, and its energy dependence interpo- 
lates between different limits. In particular, S {E) is 
quadratic in E at the lowest energies (a Fermi-liquid 
form), and is almost flat at large E. At intermediate 
energies, S {E) is roughly linear in E. The real part of 
the self-energy is linear in E at the lowest frequencies 
(— dS {E)/dEE^o ~ ^n)t s-nd is flat at high frequencies. 
If the relaxational spectrum of spin fluctuations is cut at 
some upper cutoff, the real part of the self-energy will 
start decreasing above the cutoff. 

In the superconducting state, the self-energy changes 
by two reasons. First, fermionic excitations acquire a 
gap. Second, the spectrum of collective excitations by 
itself changes as a feedback from the gap opening. The 
expression for the self-energy incorporates both effects 
and is given by 
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l-n{E')/LOsf ^{E + E' 
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where n(i?) is the polarizability bubble in the supercon- 
ducting state (a sum of the two bubbles made of normal 
and anomalous fermionic Green's functions). This polar- 
ization operator can be computed explicitly. We obtained 



n{E) 



2A^ \iE^ ) 



for \E\ < 2A 
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n (E) 



^D (^ ] for \E\ > 2A 



for \E\ < 2A 



\E\ ^2 (1 - 4^ ) for 1^1 > 2A 



(18) 



where D{x^) = {Ki{x^) - K2{x^))/x'^, and Ki{x'^) and 
K2{x^) are the elliptic integrals of flrst and second kind, 
respectively. The expression for 11 was earlier obtained 
ii)2i. 

We see that n" is flnite only at |£:| > 2A. At 
\E\ < 2A, n (E) is positive and interpolates between 
zero at E = and inflnity at \E\ = 2 A (at the low- 
est energies, !!(£') w (7r/8)-E'^/A). At some frequency 
flo, n (r^o) = ^sf, and the dynamical spin susceptibility 
Xs{E) oc 1/(1 — Il{E)/LLisf) has a pole. As a result, the 
gapless continuum of the normal state splits into two sep- 
arate entities: the gapped continuum at energies above 
2A, where 11 is non-zero, and the pole (the resonance 
peak) at an energy flo below 2 A. We see therefore that 
in the superconducting state, one-mode and continuum 
models are quite similar - both describe the interaction 
between fermions and a bosonic mode. The difference 
between the two models is in the details, and also in the 
fact that in a collective mode model, the bosonic spec- 
trum still contains a continuum above 2 A. 

The location of the pole can be straightforwardly ob- 
tained from ifTSJI . For A = 30 meV, the mode is at flo = 
40 meV, if tOsf ~ 26 meV. This last value is quite con- 
sistent with earlier estimates^l Near the pole, the spin 
susceptibility is 



X{E) « Zo 



no 



nl-E^- iS 



(19) 



where Zq ^ 0.77. Apart from the residue Zq, Eq. iflfljl 
describes the same propagator as in the one-mode model 
(see Eq. Q). 

Substituting the results for 11 into the expression for 
the fermionic self-energy, Eq. ((T6|l . we flnd S(i?) as a sum 
of two contributions. One comes from the pole and the 
other comes from the gapped continuum. The generic be- 
havior of the self-energy is similar to what we have found 
for the one-mode model. Namely, S'(i?) is linear in E at 
the lowest energies, and diverges as a square-root at ap- 
proaching — (A -I- Oq) from below. Above this threshold, 
S (E) drops to a flnite value, and decreases at even larger 
\E\. The imaginary part of the self-energy is zero below 
the threshold at — (A -I- fJo), diverges as a square-root at 



approaching the threshold from larger \E\, and eventu- 
ally recovers the normal state value at highest energies. 
At the smallest E, we found that the dominant contribu- 
tion to 5]'(i?) = — Ag^i? comes from the mode, continuum 
only accounts for about 20% percent correction. Evalu- 
ating the integrals, we found that A^^ ~ O.TAJj. This is 
similar to what we have found in the one-mode model. 
At -E = —A, we found, using the full form of the polar- 
ization bubble, A^^ ^ ~ 0.75A$^, which is again similar to 
what we have found in the one-mode model. 

For the imaginary part of the self-energy E {E) and 
i? < we found 



where 
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The first contribution is from the mode, the second is 
from the gapped continuum. At —(Sip + A) > E > — 3A 
only the mode contributes. The self-energy in this range 
is very similar to the one- mode result. Above 3 A, the 
gapped continuum also contributes to S (E), initially 
as ^/E + 3A/\og^{E + 3A) for E < -3A, and more 
strongly at larger \E\. Combining the contributions from 
the mode and from the gapped continuum, we found 
numerically that the total S is almost fiat above 3A 
at a value S sa 1.5A5JJ7o- The near-constant value of 
E is quite close to the normal state value in the one- 
mode model, S]„ = (7r/2)A^r2o, but we stress that in the 
collective mode model, this fiat behavior is obtained at 
|i?| « 3 -^ 5A. In the one-mode model, S {E) at these 
energies has a strong frequency dependence ranging be- 
tween 1.25Sl' at E = -3A to 1.04E" at £^ = -5A. 



C. Comparison and application of the two models 

Not surprisingly, one-mode and collective mode models 
give very different results for the normal state. Within 
the one-mode model, the normal state spectral function 
still shows a peak-dip-hump structure, and the renor- 
malized dispersion displays an S'-shape structure near 
E = — ilo- In the collective mode model, the imaginary 
part of the self-energy is roughly linear in E at frequen- 
cies comparable to Qq, and S (E) displays a crossover 
from a linear behavior at small frequencies to a near con- 
stant behavior at higher frequencies. From this perspec- 
tive, a combination of the measurements below and above 
Tc provides the best way to distinguish between the two 
models, particularly as we found the relation between the 



coupling constants in the normal and superconducting 
states. Several ARPES measurements near (7r,0) indi- 
cat oTiffi'? that the coupHng to the mode disappears above 
Tc thus strongly supporting the collective mode model. 

The normal state measurements may be "contami- 
nated" by thermal effects, which mask the difference be- 
tween the two models. A way to avoid thermal effects is 
to focus on low T measurements. However the two mod- 
els give very similar results for the superconducting state. 
The only qualitative difference is the gapped continuum 
which is still present in the collective mode model, but 
the continuum affects the self-energy only in a moder- 
ate extent both at vanishing E and at \E\ ~ A. The 
dominant contribution to the self-energy at these ener- 
gies comes from the resonance at SIq, which is present in 
both models. The continuum does affect the self-energy 
at \E\ ^ 3 -f- 5A, but it is difficult to measure the self- 
energy in the (tt, 0) region in this energy range since the 
bare dispersion only extends to ~ 2A for the antibonding 
band and to ~ 7A for the bonding band and therefore 
all evaluations strongly depend on the exact values of the 
bare particle dispersion. 

On the other hand, the close similarity between the 
two models in the superconducting state is good for ad- 
dressing the fundamental issue whether the data in the 
superconducting state are actually consistent with the 
collective mode model, and with estimates of A^^ or A^. 

For the analysis of the experimental data we used the 
one-mode model which is determined by A^^ and which 
simulates the coupling to the magnetic resonance modei^ 
The bare one-mode model is extended by adding to the 
self-energy a Fermi-liquid-like term which is shifted by 
3A to higher binding energy. This term approximates 
the gapped continuum and is determined by the cou- 
pling constant A^^,. We then compare the two coupling 
constants A^^, and A^ with theoretical estimates for the 
coupling constant A^c in the collective mode modeli^i*^ 
We obtain reasonable agreement between experiment and 
theory which indicates that the dressing of the charge car- 
riers in the (7r,0) region is related to a coupHng of spin 
excitations. We also show that the magnitude S"(— oo) 
evaluated from the incoherent spectral weight is consis- 
tent with the value of A^ derived from the dispersion near 
A. This means that the spectral function for the coher- 
ent and the incoherent states can be described by one 
self-energy function indicating a common linear dressing 
for the both states. 



III. 



EXPERIMENTAL SETUP 



The ARPES experiments were carried out at the 
BESSY synchrotron radiation facility using the U125/1- 
PGM beam line and a SCIENTA SESIOO analyser. Spec- 
tra were taken with various photon energies ranging from 
17 to 65 eV. The total energy resolution ranged from 8 
meV (FWHM) at photon energies hj/=17-25 eV to 22.5 
meV at hj^=65 eV. The momentum resolution was set to 




0.05 



WAVE VECTOR (A"' 



-D.I 1 
WAVE VECTOR (A'l 



Figure 4: ARPES intensity plots as a function of energy and 
wave vectors along the (tt, n) — (tt, — tt) direction of the opti- 
mally doped Pb-Bi2212 superconductor taken at T = 30 K. 
Zero corresponds to the {it, 0) point, (a) bonding band, (b) 
antibonding band, (c) and (d): calculated spectral function 
using a model self-energy function for the region around the 
(tt, 0) point. The red dashed line represents the bare-particle 
bandstructure 



0.01 A^^ parallel to the (tt, 0) — (tt, tt) direction and 0.02 
A~^ parallel to the T — (tt, 0) direction. Here we focus on 
spectra taken with photon energies of 38 eV and 50 (or 
55 eV) to discriminate between bonding and antibonding 
bands. The polarization of the radiation was along the 
r — (tt, 0) direction. Measurements have been performed 
on (1x5) superstructure-free, optimally doped BiPb2212 
single crystal with a Tc= 89 K. Since data for the nor- 
mal state have already been published^ we only show 
data measured at T = 30 K. 



IV. EXPERIMENTAL RESULTS 

In Fig. 21 we show typical ARPES data for wave vec- 
tors close to the (tt, tt) — (tt, — tt) line, centered around 
the (tt, 0) point. As has been shown previousljiStiSiii and 
supported by theoretical calculation a^"'^'^ , the data taken 
with hw=38eV due to matrix element effects represent 
mainly the bonding band with some contributions from 
the antibonding band. The data taken at h'L'=50 eV (see 
Fig. ^b)) have almost pure antibonding character. In 
order to obtain the spectral weight of the pure bond- 
ing band (see Fig. El a)), a fraction of the 50 eV data 
has been subtracted from the 38 eV data. We have also 
added in Fig. 31 the bare-particle dispersion which was 
obtained from a self-consistent evaluation of the data 
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Figure 5: ARPES intensity plot for k- values near the 
{lAn, n) - {lAn, -tt) line of the optimally doped Pb-Bi2212 
superconductor taken at T = 30 K. Zero corresponds to the 
(1.47r, 0) point. The data were taken with a photon energy hv 
= 50 eV in order to maximize the intensity of the antibonding 
band. 



at the nodal poinfe^S, an evaluation of the anisotropic 
plasmon dispersion^LS and from LDA bandstructure 
calculationsSS. When comparing this bare-particle dis- 
persion with the very broad distribution of the bond- 
ing band at high energies, one realizes a renormaliza- 
tion of the occupied bandwidth by a factor of about 
1.7 corresponding to A^^ « 0.7. This value is not far 
from that derived for the bandwidth renormalization 
above Tc, where no additional renormalization effects 
at lower energies have been detected^ and from that at 
the nodal poinUSiSi. Probably a large fraction of this 
bandwidth renormalization stems from a coupling of the 
charge carriers to the above mentioned continuum of spin 
fiuctuationsSi. 

In Fig. we show an ARPES intensity distribution 
near kp of the antibonding band, close to the (1.47r, n) — 
(1.47r, -tt) line, of Pb-Bi2212 measured at 30 K with a 
photon energy hv = 50 eV. At this place in the second 
Brillouin zone, the bare particle dispersion of the an- 
tibonding band reaches well below Eb = 70 meV and 
therefore contrary to Fig. 4 (b) the branching into two 
dispersive arms can be clearly realized. These data to- 
gether with the data of Fig. when compared with the 
model calculations shown in Fig. clearly reveal that 
the dominant effect of the renormalization, besides the 
bandwidth renormalization mentioned above, is due to a 
couphng to a bosonic mode leading to a branching energy 
of -- 70 meV. 

In order to obtain more quantitative information on 
the parameters which determine the self-energy function 
leading to this renormalization we have performed vari- 
ous cuts of the spectral weight shown in Fig. 01(a) which 
are presented in Fig. El A constant-k scan for k=kF is 
depicted in Fig. EJa) showing the typical peak-dip-hump 
structure presented for data taken at the (tt, 0) point in 
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Figure 6: (a) Constant-k cut of the data shown in Fig at 
kp; (b) constant-k cut of the data at about one third of kp 
(starting at the (tt, 0) point); (c) constant-E cut at E = -100 
meV, The solid lines represent fits to the data (see text) 



previous studiesiiSii. From the peak energy one can de- 
rive the superconducting gap energy A = 30(4) meV. 
In the previous HteraturaS^, from the dip energy at ^ 
70 meV the branching energy Eb was derived. Another 
constant-k scan at 1/3 kp (starting from the (tt, 0) point) 
is shown in Fig.EJb). At this k- value the intensity of the 
coherent peak is strongly reduced and in the framework 
of the one-mode model, mainly the threshold of the in- 
coherent states (the hump) is observed. From a fit to 
these data, taking into account a small intensity of the 
coherent line and a threshold of the incoherent states, the 
threshold energy could be determined which in the one- 
band model yields the branching energy, Eb = A + Qq 
= 70(5) ( meV. This together with the gap energy A= 30 
meV yields a mode energy of 40 meV. 

In Fig. ^c) we show a constant-energy cut aX E = - 
100 meV of the data presented in Fig, As discussed 
in Sect. ^ from the fit of those cuts with a Lorentzian 
one can obtain from the width of the Lorentzian a value 
of the imaginary part of the self-energy at the selected 
energy. For the energy below the mode energy, this value 
E"(— oo) is a measure of the coupling to the mode. The 
actual situation is more complicated since near the antin- 
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Figure 7: Typical values for S"(— oo) as derived from a fit of 
constant-E cuts (see Fig, EJ of data shown in Fig, Ela) and 
(b) using Eq, (4) The solid line is a guide to the eyes, (a) 
bonding band, (b) antibonding band. 



odal point, the bandstructure is far from being Hnear at 
this energy range. We have fitted the data shown in 
Fig. ISIc) by Eq. Q with the self-energy given by jSJ, 
the bare dispersion extracted from earlier work^^ A = 
30 meV, Qq = 40 meV, and E"(— oo) used as a param- 
eter. Moreover as described above a Fermi-liquid like 
term was added to the self-energy to approximate the in- 
fluence of the gapped continuum. The imaginary part of 
this term is given by a{E — 3A)^ for \E\ > 3A and zero 
for \E\ < 3 A where the magnitude of a is determined by 
A{g. The value 3 A can be easily understood by looking 
at Fig. n since in the superconducting state the states 
available for a fermion decay have minimum energy of 
3 A, and hence below 3 A, the correction to E from the 
continuum is zero. Using this self-energy function from 
an extended one-mode model the flt yielded typical val- 
ues for the parameter i;"(— cx)) as shown in Fig. For 
\E\ > Eb = 70 meV the values should be constant. The 
finite slope detected in the analysis may be related to 
errors in the bare particle dispersion, to the assumption 
of a constant density of states during the definition of 
the self-energy function, or to in the assumed A^^. For 
\E\ < Eb the results from the fits are determined by the 
fiat dispersion of the coherent states and therefore, due 
to the finite energy resolution, large values are obtained 
in this energy range (not shown). From evaluations of 
such data taken on several samples we derive a value 
S"(-oo) = 130(30) meV. The large error for this value 
stems from various measurements on samples with slight 
mismatches in there orientation which leads to different 
bare bandstructures as compared to the assumed one. 

Another important information comes from the dis- 
persion of the coherent spectral weight between the gap 
energy —A and the branching energy — (A -f Qq). Origi- 
nalljiSii, the data were fitted using Eq. (fT3ll . As pointed 
out in Sect, m this is a good approximation for conven- 
tional superconductors, where the mode energy is much 
higher than the gap, but not for the high-Tc supercon- 
ductors. Here the energy of the mode is comparable to 
the gap energy and therefore the renormalization func- 
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Table II: Parameters determining the self-energy function 
of Pb-Bi2212 near {{tt,0) below T^. A-'', A'', A*, coupling 
constants from bandwidth renormalization, from coupling to 
a bosonic mode, and total coupling constant, respectively. 
E"(— oo) is determining the imaginary part of the self-energy 
function at high binding energies, sc: parameters derived in 
the superconducting state, n: parameters calculated for the 
normal state by setting A to zero. The data are compared 
with parameter derived for a Mo(llO) surface stated and for 



system 



A-'' 



A" 



A* 



E"(-oo)(meV) 



BiPb2212 sc 


0.7(3) 


2.0(4) 


2.7(5) 


130(30) 


BiPb2212 n 


1.3(3) 


2.7(4) 


- 


130(30) 


Mo(llO) 


- 


0.4 


- 


15 


Pb n 


- 


1.6 


- 


- 



tion from which the A^^ -values are derived, depends on 
energy and also on A/^o- In this communication we 
have fitted the data using the full Eqs. Q) and ltT2)l . 
For the complex self-energy function, we have used the 
same parameters as for the fit of the constant-energy cuts 
shown in Fig. |Sl(c). We then obtain a complex renor- 
malization function Z(E) from which we derive A- values 
not at energies between A and A -I- J7 as in the previ- 
ous studji^ but at zero energy, i.e., A^j, = Z(0) — 1. In 
the fit we have used the above given values for A and 
r^o and we have chosen A^^ as a parameter. From the 
fit, we extracted A^^ = 2.0(4). This yields a total Xlc= 
2.7(5) from Z(0) in the superconducting state composed 
of a bandwidth renormalization part A^^, = 0.7(3) and 
a bosonic part A^^, = 2.0(4). The errors are not due to 
statistics but result from different measurements on dif- 
ferent samples. Using the relation A^^ = 0.74A^ (see Sec. 
ITTjl we obtain Afj = 2.7. To obtain A^ in the normal 
state we used the same ImT, from fermion-fermion inter- 
action as in the superconducting state, but set A = 0. 
This yields A^ ^ 1.3 and a total coupling constant for 
the normal state A^ = 4.0(5). We cohected the values 
of the coupling constants and S]"(— oo) in Table II. Note 
that the normal state values are not derived from mea- 
surements in the normal state. Rather they were derived 
from data taken in the superconducting state and setting 
A to zero in the renormalization function. This means 
that A^ in Table II is a fictitious normal state coupling 
constant because the bosonic mode does not exist in the 
normal state. It is only presented for the comparison of 
the coupling constant in HTSCs with those of conven- 
tional metals and superconductors. 



V. DISCUSSION 

The coupling constant to a bosonic mode was de- 
rived from the dispersion of the coherent states between 
-30 and -70 meV while S"(— oo) was derived from a 
constant-E cut of the incoherent states below -70 meV. 
If the data can be described by one self-energy function 



which essentially results from a coupling to one bosonic 
mode, the A^-value calculated from !]"(— oo) within this 
model should agree with the A^-value extracted from 
the dispersion of the coherent states. Using the rela- 
tion A^ = -2E"(-oo)/(7rrJo) (see Sect.|n|| and the value 
S"(-oo) = 130(30) meV one obtains Aj; = 2.1(5) which is 
in reasonable agreement with A^ ~ 2.7(4) derived from 
the dispersion near —A. This supports the idea that 
the coherent and the incoherent spectral weight of the 
spectral function in the superconducting state can be de- 
scribed by a single self-energy function, which is essen- 
tially determined by the coupling to one bosonic mode at 
40 meV. This view is also supported by the fact that tak- 
ing this self-energy function and calculating the spectral 
function using Eq. JSJ we obtain a reasonable agreement 
with the experimental ARPES data, both for the bonding 
and the antibonding band (see Fig. 0)1 . The differences 
in the intensities of the bonding band near (7r,0) may 
be explained by matrix element effects. We emphasize 
that in the superconducting state both the real and the 
imaginary part of the self-energy function indicate a very 
strong couphng to a bosonic mode. 

Recently, there has been some evidence from ARPES 
measurements^iiS^ that in the undoped cuprates there 
is a very large electron-phonon coupling leading to a 
strong polaronic renormalization connected with a negli- 
gible spectral weight for the coherent states and a high 
spectral weight for a multiphononic line at higher ener- 
gies. Furthermore there are theories of the pairing in 
high-Tc superconductors which are based on the forma- 
tion of polarons and bipolarons^. The big question is 
whether this strong electron-phonon coupling survives 
for the high-Tc superconductors or whether it will be 
screened by the charge carriers and at what dopant con- 
centration the adiabatic approximation is valid, where 
the Fermi energy is much larger than the mode energy. 
The data shown in Fig. which were actually taken 
down to an energy of -400 meV, show no indication of 
a polaronic line at lower energy. Moreover, the data in 
Fig. and its evaluation in terms of a coupling to a sin- 
gle bosonic mode gives no room for multi-bosonic polaron 
excitations for optimally doped samples. 

The analysis given above clearly identifies below Tc a 
very strong couphng to a bosonic mode. In Table II we 
have listed other couphng constants to bosonic (phonon) 
modes detected by ARPES. Compared to the surface 
state of the Mo(llO) surface coupled to a phonon mode, 
A^ in OPBiPb2212 is a factor of 6 larger. Almost the 
same factor 8 is obtained for S"(— oo). Compared to 
the strong coupling superconductor Pb the coupling con- 
stant A^ for OP BiPb2212 is a factor of 1.6 larger. This 
indicates that we really have a very strong couphng to a 
bosonic mode. This coupling is even enhanced at lower 
dopant concentration^. 

In the following we compare the experimental coupling 
constants A^ and A^^ with those derived in the cohective 
mode model for the gapped continuum and for the single 
mode, respectively. In the collective mode model, in the 
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superconducting state about 20 % of the total coupling 
constant comes from the gapped continuum (see Sec. II 
B). The corresponding experimental value A^/A*^, = 0.26 
(see Table II) is in remarkable agreement with the the- 
oretical value. This also holds for the absolute values of 
the coupHng constant. In previous work^LSi the normal 
state coupHng constant for the collective mode model was 
estimated to be A^ = 2 -^ 3 from fits of theoretical values 
for ujsf, r^o, and Tc. This transforms into A^^, = 1.5-f-2.25 
for the superconducting state which is not far from the 
experimental value A*^ = 2.7(5). Thus the agreement 
of the relative and absolute experimental values of the 
coupling constants with those derived from the collective 
mode model is a strong indication that the dressing of the 
charge carriers in the (7r,0) region (i.e. the region were 
the superconducting order parameter has its maximum) 
is predominantly determined by a coupling to spin exci- 
tations, in particular to the magnetic resonance mode. 

This interpretation is supported by several other 
ARPES results. The strong temperature dependence of 
the coupling to the modeS*^ is difficult to understand in 
terms of electron phonon coupling. Our model calcula- 
tions also show that the data above Tc cannot be de- 
scribed by a thermally broadened phonon line. Also the 
strong dopant dependenceSii is difficult to be explained 
in terms of electron phonon coupling. Furthermore there 
is a large coupling at (tt, 0) and a much smaller coupling 
to the mode at the nodal point, indicating a coupling of 
states which are separated by a wavevector (tt, tt) typi- 
cal of an antiferromagnetic susceptibility. Moreover the 
energy of the bosonic mode detected in ARPES is close 
to the energy of the magnetic resonance mode detected 
in inelastic neutron scattering. Finally we mention re- 
cent ARPES measurements on the parity of the coupHng 
between bonding and antibonding band^ and the "mag- 
netic isotope effect", i.e., the strong changes of the dress- 
ing of the charge carriers upon substitution of Cu by 



Zn4^i^, which both support the magnetic scenario. 



VI. 



CONCLUSION 



In this contribution we have analyzed the spectral 
function of optimally doped BiPb2212 near the antin- 
odal points measured by ARPES. Compared to previous 
studies, we have not analyzed just one constant-k cut 
or just the dispersion of the coherent state but the 
entire spectral function including the coherent and the 
incoherent spectral weight. In this context we have used 
expressions which not only can be used in the case of 
normal superconductors but also for HTSCs where the 
mode energy is not much larger than the superconduct- 
ing gap. It was possible to describe the spectral function 
using a single parameterized self-energy function. By 
comparison of the experimental data with theoretical 
models, we conclude that the main contribution to the 
self-energy is a very strong coupling to the magnetic 
resonance mode. At higher energies (and above Tc) it 
was necessary to take into account a bandwidth renor- 
malization by a factor of two due to interaction with a 
gapped (ungapped) continuum of spin excitations. There 
is no evidence for multi-bosonic polaron excitations for 
this dopant concentration. 
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